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1 Introduction: periodic media nomenclature 

1. The space domain is defined by a basis, (a i, a,2, as), where any vector can be written as 

f' = r + R = f + oiai + a 2 a 2 + 0303 , (1) 
where R is the translation vector, with a±, a.i, 03 integers. 

2. The spectral domain is defined by a basis, (61, 62, 63), and similarly, the translational vector 
is written as 

G = ^61 + (3 2 b2 + Psh , (2) 

where (3\,P2,Ps are integers. 

3. The two basis are linked since the functions (fields, permittivity) are periodic. For example, 
if we write the permittivity: 

Fourier expansion: e(r) = e(G) e tG ' r where e(G) = — — jjj dr 3 e(r) e~ tG ' r . (3) 

_ cell 

Periodicity: e(r + R) = ^ e(G) e iG < f+R ^ 

G 

= ^ e(G) e^^ R = e(f) (4) 
G 

so that R = 1 and 

G-R = 2mir where m €{..., -1,0, 1,2, ...} . (5) 

We can see that condition (5) is immediately verified if we impose: 

bj • Hi = 2wSij . (6) 
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Section 1. Introduction: periodic media nomenclature 



4. Bloch-Floquet theorem: 

Since EM fields are periodic, we can write them as a propagating function times a function 
with the same periodicity as the medium: 



= e^ f C£(r) where (^f + R) = ^(f ) , 

and where £ can represent either the electric or magnetic fields, E or H. 
Since ((f) is periodic, we can Fourier expand it: 



(7) 



so that we shall write: 



5. Wave equation in source-free region: 



G 



G 



:(k+G)-f 



^(r) = E^ e4+G) ' f 



(8) 



(9a) 
(9b) 



From Maxwell's equation, we can easily obtain the following wave equations in source-free 
regions (with e = e(r)): 



V x 



V x V x E{f) = ( ^ /fx r e r (r) £(f) 
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e r (r) 



— V x H(f) 



- ) Mr H(r) , 



(10a) 
(10b) 



To make these equations more symmetrical, we shall work with l/e r (r) instead of e r (r) 
directly, so that we define 



K r (r) = 



1 



e r (r) 



J>(G). 



tG-r 



(11) 



G 



The wave equations are rewritten as: 



K r (r)V x V x E(f) 



V x 



K r (r)V x i?(r) 



UJ 



UJ 



Hr E(r) , 



(12a) 
(12b) 
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2 Treatment of the E field 

2.1 Method 1: direct expansion of the permittivity 

We want to write Eq. (10a) with the decomposition of Eq. (9a). First, let us compute the first 
curl (taking G' as the variable for the expansion): 



e G' e 



Vxfi fc (f) = ^Vx 
G' 

Taking the curl one more time gives 



(k+G )-f 



iJ2(k + G')xe^ e ^ k+G > r 
G' 



(13) 



VxVx£ I (f) = -3H(/) 



G 



(k + G) x e^/ 



i(k+G )-r 



Upon using Eq. (3) but changing the index G into G , we write 



e r (f)E(f) = ^2Y1 e ~^") 5 G' el{k+G +G ^ ■ 
G' G" 



By changing the variables G = G' + G" ': 



(14) 



(15) 



e r (f)E(f) = £ £ ~ e r(G ~ G') eg e*^* . 
G G' 

The wave equation (see Eq. (10a)) can therefore be rewritten as: 



(16) 



g' 



(k + G ) x e 



G 



e KMyf = (^Y^J2J2e r (G-G , )ege^ k+G > r . (17) 
G G' 



We can simplify by exp (ik • r) and multiply by exp (— iG • r) to get: 



G' 



(k + G)xe 



G 



e« G '- G ">r= fc) ^E^G-G')-e^e^ G - G '^. (18) 



G G 



If we integrate this equation over the entire space, we can pull all the terms out of the 
7=»' t=<" — t=? 7=7" — 

integral, except ~ • > ' r on the left-hand side and. - > r on the right-hand side. Yet, 

we have 



v 



d ^3 e^ (G'-G").r = _l_ ( 5 ( G-G ,, ), 



(19) 



so that Eq. (18) becomes (upon substituting G by G since these are dummy variables): 



-(k + G) x 



(k + G)xe- a = (^j n r ^r(G - G')e^ , VG. 

G 



(20) 



2.2 Method 2: expansion of the inverse of the permittivity 



2.2 Method 2: expansion of the inverse of the permittivity 

Instead of working with Eq. (10a), we can also use Eq. (12a), which would need the expansion 
of Eq. (11). 

Applying the same method (and transforming the index of Eq. (11) from G to G ), we get: 



-EE^ G ")( fc+G ') 



g" g' 



(k + G ) x e 



G 



e i(k + G' + G"yr = (^Y/x r £e^e^ +G >. (21) 

G' 



which, upon substituting G = G + G , simplifying by exp(ik-f), multiplying by 
exp(— iG" -f), integrating over the whole space, using Eq. (19) and finally substituting G" 
by G, becomes: 



J2^r(G-G')(k + G') 
G' 



(k + G)xe 



G 



LO 



G 



VG. 



(22) 



3 Treatment of the H field 

The H field is treated in an exactly similar way to eventually obtain very similar equations. 

However, these equations can still be pushed further by using the fact that V • H-^(r) = 0. 
Upon using this equality, we see from Eq. (9b) that (using G' for the expansion of the field): 



(k + G ) ■ = . 
G 

We can therefore define three vectors (ei,e2,es) such that 

k + G' = \k + G'l e 3 , 
ei • e 3 = e 2 • e 3 = , 

and (ei,e2,es) for an orthonormal tryad. In that case, we can decompose 

hp = h X pe\ + h 2 pe 2 = ^ h X p e x . 



A=l,2 



(23) 

(24a) 
(24b) 

(25) 



We need now to introduce this expression into Eq. (12b). First, we compute 



G 



G 



(k + G)xe x 



i(k+G )-r 



so that 



G" G' A 



(HG)xe A 



3 i(/c+G +G )-r 



= * EEE^G' «;r(G - G' 



G G' A 
Taking the next curl, we write: 



(HG)xe A 



i(k+G)-r 



(26) 



(27) 
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Vx 



J G G' x 

so that the wave equation (see Eq. (12b)) becomes: 



(k+G)x[(k+G)xe x ] 



i(k+G)-f 



(28) 



G G 



(k+G)x[(k+G)xe x ] 



G 



(29) 

Always by the same token (multiplying by the proper functions and integrating over whole 
space), we write: 



£E*v Mg-g') 



G 



{k + G)x[{k + G)x e x ] 



C ' X" 



We can further simplify this expression by dot-multiplying the equation by ey and noting 
that (using C ■ (A x B) = B ■ (G x A)) 



(k + G)x[(k + G)x e x ] 



(k + G ) x e x 



(k + G)x e w 



(31) 



Therefore, dot-multiplying Eq. (30) by ey, we get the final result: 



EE 



G 



{k + G) xex 



{k + G)x ey 



K r (G - G) hx^i 



r f Vrhy-Q- (32) 



Upon exchanging G and G' (transformations: G — > G" , Cf — > G, G" — > G'), we obtain 



EE 

G A 



(fc + G) x e A 



(HG)xe A - 



K r (G - G) /i A £ = ( — ) Hrhx'rii ■ (33) 



G 



which is the relation given in [Joannopoulos et al, 1995, p. 129]. Upon using the same 
notation, we rewrite Eq. (33) as: 



aG 



k h ( UJ \ 

(xGuxGy \xG) = { ~ c ) » r \xGy ' 



where 
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k_ 
(aG),(aG) 



, = «r(G - G) 



(fc + G) x e A 



(HG)xey 



(34a) 



(34b) 
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3.1 Matrix form 



3.1 Matrix form 

We can cast Eq. (33) in matrix form. First, we rewrite the kernel of the operator of Eq. (34b) 
as: 



(k + G) x e A 



(k + G ) x ey 



{k + G) 



{k + G 



[e 3 x e A j • e 3 x e x > . 



Remembering that e 3 x e\ = e2 and e 3 x e2 = — e±, we can write: 



[e 3 x e A ] • [e 3 x e A '] 



e2 • e2 — e2 • ei 
-ei ■ €.2 ei • ei 



so that we write the operator as: 



(35) 



(36) 



(aG),(aG) 



,=«r(G -G) 

= R r (G' - G) 



(k + G)x e x 



{k + G) x e\> 



{k + G) 



{k + G 



e2 ■ e2 — e2 • ei 
-ei-e 2 ei-ei 



(37) 



used in Eq. (34a) 



